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Abstract 

We show how the two-body potential may be uniquely determined from n-body 
spectra in cases where the hypercentral approximation is valid. We illustrate this 
by considering an harmonic oscillator potential which has been altered by changing 
the energy or normalisation constant of the ground state of the n-body system and 
finding how this modifies the two-body potential. It is shown that with increasing 
number of particles the spectrum must be known more precisely to obtain the two- 
body potential to the same degree of accuracy. 
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1 Introduction 



The standard inverse scattering theory is formulated for two-body systems Its aim is 
the determination of the interaction from the spectral information pertaining to bound 
and scattering states of the two particles. In some cases there is no or insufficient spectral 
information available on the two-body system, while there is access to the spectral data of 
a corresponding many-body system. An example of this is the baryon spectrum whereas 
no two-quark system is known to exist. Thus it is equally interesting to ask whether it is 
possible to extract the two-body potential from the bound and scattering states of a many- 
body system. To our knowledge this general question has never been considered. There 
is a schematic treatment of the three-body problem in one space dimension by Zakhariev 
0, but this is unsuitable for practical applications. With respect to applications, only 
recently M has a feasible method for the inversion of baryon spectra been proposed and 
have calculations of the associated quark-quark potentials been reported. 

In the following, spectral information given in terms of many-body bound states is 
considered. For many-body bound states there are approximations which bear a formal 
resemblance to a two-body Schrodinger equation. This property suggests the possibility 
of applying standard two-body inversion techniques within these approaches. Examples of 
these include the Hartree-Fock method in which the interaction of the bodies is described 
by a mean field, and the hypercentral approximation (HCA) in which only the hyperradial 
part of the interaction in the hyperspherical harmonic expansion method || (HHEM) is 
retained. It is on the latter scheme, described in more detail in the following section, that 
we will focus in this article. The hypercentral potential used in this method is determined 
uniquely by the underlying two-body interaction. The main progress achieved in this 
study is an inversion of this relation for a general n-body system consisting of n identical 
particles. This inversion procedure is formulated for arbitrary n and can be applied 
as long as the HCA gives a proper description of the many-body system. In principle 
the method is an extension of the inversion procedure of baryon spectra, investigated in 
previous papers @, ]5[. In this extended scheme we show how the two-body interaction 
may be obtained in general from n-body spectral information, and in particular use the 
method to obtain the two-body interaction from three-, four-, five- and six-body spectral 
data. 

We illustrate this scheme by applying our inversion method to spectra which coincide 
with that of the harmonic oscillator in both energy levels, and normalisation constants 
of the bound states in hyperspherical space, bar the ground state. This ground state is 
either shifted in energy relative to the harmonic oscillator ground state or has a different 
bound state normalisation constant than that of the harmonic oscillator ground state. 
Primarily we focus on the changes of the two-body potential with increasing number of 
particles n. Our example clearly demonstrates that the accuracy of the determination 
of the two-body potential from the n-body spectrum deteriorates with increasing n, if 
the accuracy to which the bound states are given remains fixed. It implies a rapid loss 
of information on the underlying two-particle potential with increasing size of the bound 
system of particles, despite the corresponding increase in the number of states present 
in their spectrum. This can also be turned around and viewed as an indication that the 
n-body spectrum becomes less sensitive to the details of the two-body interaction with 
increasing number of particles. 

The details of the inversion in the HCA are described in the following section, and 
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the illustration and the results are to be found in section 3. We conclude in section 4, 
discussing what may be learnt from the example we have used. 

2 The hypercentral approximation and inversion for- 
malism 

We begin this section with a description of the HCA and in particular show the relation 
of the hypercentral potential to Wigner-type two-body interactions. We then show how 
this relation may be inverted, thereby facilitating an inversion of n-body spectral data to 
obtain the underlying two-body interaction. 

2.1 The hypercentral approximation 

An n-body state is described after removing the centre-of-mass motion by 3(n-l) internal 
co-ordinates. The corresponding n-particle Schrodinger equation can generally only be 
solved using various approximations. A very effective procedure is the expansion of the 
wave functions and potentials in terms of hyperspherical harmonics. This HHEM is 
particularly useful for the evaluation of many-body bound states of both bosons and 
fermions. 

Specifically for fairly soft two-body interactions, where correlations between the parti- 
cles are negligible, the lowest order of the HHEM, the so-called hypercentral approxima- 
tion (HCA), provides an excellent description of the many-body bound states. It is exact 
in the case of the n-body harmonic oscillator, where it is well know that the n-body bound 
states may be obtained as linear combinations of single particle oscillator eigenfunctions. 
A prominent example of the applicability of the HCA is the calculation of baryon spectra 
H in the framework of the non-relativistic quark model, in which the forces operating 
between the constituent quarks are soft and confining. Even for a Coulomb interaction or 
certain nucleon-nucleon interactions, the HCA alone may account for up to 90 % of the 
binding energy of an n-body bound state of the system (see |2[] and references therein) . 

In the HCA the n-body Schrodinger equation reduces to an effective Schrodinger 
equation in the hyperradius p 

o n 

P 2 = -E(r,-r,) 2 , (1) 

where denotes the co-ordinate of the z-th particle. The hypercentral potential V^ c (p) 
which accounts for the interactions of the n bodies is the rotationally invariant part of 
the sum of two-body central potentials C/(|rj — tj\) acting between each pair of particles 
i and j in the 3(n — l)-dimenional space. It is determined by the two-body potential U 
through the linear relation M, ff\ 

v h n c ( P ) = E^ n (p) 

3=0 

= ^S^dz{l-z)^ L ^{l + zy^U{ps^l (2) 

where D has been used to denote 3(n — 1), and l m is the maximum value of the orbital 
angular momentum in the outermost shell consistent with a hyperspherical harmonic of 
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degree L m . By allowing the minimal degree L m of hyperspherical harmonics present in 
the expansion of the n-body wavefunction to be greater than zero, the Pauli Principle may 
be implemented for bound states of fermions. The normalisation constant h is defined by 



21 n 



-i 



h = Y / a J dz(l - z) V+^-i(l + z )3+\. (3) 

The Oj are constants determined by the shell structure of the bound state. If L m = 0, as 
is the case for all n-boson bound states, the relation simplifies to 



r(f) /■+! JM ^d=6 „ , _ A i rr/ _ /l + z, 

0F2°/ 2 - 2 r( 



KciP) = T^HTi^TD^ / dz(l - + z )*U(pJ-^). (4) 



Thus V£ c {p) represents some weighted average of the two-body potential over the range 
(0, p), reflecting the range of the possible separations of pairs of particles in the n-body 
state of hyperradius p. The total n-body interaction due to the n( - ra ~ 1 ' > pairs of particles 
interacting via the two-body interaction U is then accounted for by n ^ n ~ 1 ^ V£ c (p). For n 
particles of equal mass m, a bound state with energy E is determined in the HCA by the 
solution of 

in 2 



m 



(P_ £(£ + 1) 

dp 2 p 



2 



^K^KciP) - E\ u{p) = (5) 



where the so-called grand orbital momentum £ depends on the orbital angular momentum 
I and the dimension of the space, 

£ = I + L m + (6) 

For all two-body potentials which are less singular than r~ 3 at the origin, the hypercentral 
potential is well-defined. This implies that the hypercentral potential is also no more 
singular than this, as the leading order behaviour of the hypercentral potential is that of 
the two-body potential near the orgin. This condition does not in reality restrict the use of 
the HCA, as an attractive potential must be less singular than r~ 2 if the spectrum is to be 
bounded from below, while the presence of a repulsive singularity in the potential would 
signal the possible inappropriateness of the HCA, as correlations between the particles 
could be important. 

In physical systems where the HCA provides a good description of the spectrum, n- 
body spectral data could be used via standard inversion techniques to obtain the effective 
interaction 

£{£ + !) m n(n-l) 

p2 + h 2 2 VhAP) 17 j 

appearing in equation (Q). The corresponding two-body potential may then be deduced 
from the hypercentral potential via a further inversion step. 

In the following we present such a procedure. As the expression is more complex if 
L m > 0, we first derive the relation for an n-boson hypercentral potential, and then use 
these results to generalise to the n-fermion case. 
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2.2 The inversion of the n-boson hypercentral potential 



Here we consider the inversion of the hypercentral potential in the simplest case, where 
L m = 0. Using the new variables 

•1 + r 



x 



p 2 and y = p 2 



we rewrite the relation ([J) in a form which we can readily invert to obtain the two-body 
potential U in terms of the hypercentral potential V£ c , 



D-3 



2 r(f) 



dy^fy[x-y) 2 U(y/y). 



(9) 



We consider two distinct cases, depending on whether the system consists of an odd or 
even number of particles. 

1. If n is even, we may simply differentiate the above expression to obtain the two- 
body potential in terms of the hypercentral potential. Setting D — 5 = 2k, where 



2, 5, 8, we find 



U(Vy) 



71 



1 



1 d k+1 



2 r(* + §) Vydy k+1 



y k+ -*v h n c (Vy) 



(10) 



2. If n is odd, we write D — 5 = 2k + 1, where k = 0, 3, 6, and on differentiating 
expression (H) we obtain 



k+l 



dx k+l 



7T 



dy 



VyU(Vy) 



This is in the form of the Abel integral equation J' 

/(*) 



9{s) = / dt- 
Jo 

which may be uniquely inverted to obtain 

dg 



/(*) = - fds- 

TT JO I 



1 



1<?(0) 



(12) 



(13) 



ds y/t — S 71 \fl 

The Abel transform is a special case of the Riemann-Liouville fractional integral 
Thus we can invert ([TT| ) to obtain 



2T(k + 3)^y 
y d k+2 



o dx 



fc+2 



dx 1 ,. 
H hm 



fc+2-r/n 



VT C (^) 



The derivation of the two-body potential from the three-body hypercentral potential 
presented previously [HI [5] is an example of the above transformation with k = 0. 

While these relations above show that it is in principle possible to extract the two-body 
interaction from a hypercentral potential describing a system of arbitrary n, it is clear 
from the form of these relations that in practice the larger the number of particles, the 
more difficult it is to perform this inversion numerically as higher order derivatives enter 
the relation. We will demonstrate this quantitatively in the following section, following 
the generalisation to the n-fermion case. 
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2.3 The inversion of the n-fermion hypercentral potential 

Using the variables x and y as defined in®, we may re-express each ^"(p) appearing in 
(D as 



n- / 9 \ ~~ 2 — + 2 /-a; n _e i 



(15) 



Then VJ 1 may be found in terms of the two-body potential as described in the preceeding 
section. 

If there are an even number of fermions, setting D — 5 = 2k, where k = 2,5,8,..., we 
find 

" r .k+L„ ' 3 - 



A j y 1+ *U(Jy) 



(jyk+Lm-j + l 



V 



vnvv) , j = o,i,. ..2i r , 



where we have defined the constant Aj as 



= 2fc+Lm+ | ^_ u± 

h 



(16) 



(17) 



Taking the j th derivative with respect to y of each of the above expressions and then 
summing over j we find the relationship from which U may be determined from a known 
hypercentral potential: 



21 ,i 



d' 



d 



k+L m +l 



h-\-L rt 



} v h n c (Vv) 



dyk+L m +l 

For an odd number of fermions, writing D — 5 = 2k + 1, where k = 0,3,6,..., we have 



rfk+Lm-j + l fX qii + hjJi fu\ 

r k+L - +2 vnv^)] = b 3 / d y y zK y) 

3 J Jo yjy-x 



X 



(19) 



d x k+L m -j+l 

where the constant B d is defined as 

r(fc + f + L m -j)2^f^ a 3 

B i = 7? T (20) 

If Vfr c (x) is no more singular than x~ 2 , then use of the Abel transform once again leads 
to the expressions 



B y*U(jy) 



7T Jo dy 



k+Lm+2 



X 



k+Lm+2-r/n 



KT(v^) 



dx 



lim 



d k+ 



Lm + l 



and 



7Ty/y x ~*° dx k+Lm+l 

1 fV (jl k + L m-j + 2 



x k+L ™ +2 V n (^/x) 



(21) 



7T 



fjk+Lm-j+2 



X 



k+L m +2\rn 



dx 



j = l,2,...2/ r 



(22) 



Taking derivatives with respect to y and performing an integration by parts we find 



B, 



d j 
'dyl 



y J+ *u(Vy) 



+ 



1 rv d k+Lm+2 r k+L +2jrn/ dx 

1 1 (jk+Lm+l 



7T ^/y — X C?X fc+Lm+1 



x fc+Lm+2 V-"(v^)l U =0 ,j = 1,2,. ..2/^23) 



Summing these expressions over j we obtain a differential equation 



2l m Jj 

U d y J 



y j+1 *U(^y) 



IT JO 



d k+ 



L m -\-2 



+ 



(lnrk-]-L r]rl -\-2 



X 



k+L m +2 



+1 



x 



k+L m +2 



lie. 



{Vx)\ 



dx 



x=0 



(24) 



as for the even number of fermions. 

Thus in order to obtain the underlying two-body potential from an n-fermion hyper- 
central potential, a differential equation of order 2l m must be solved. Note that if L m = 0, 
as we have not considered non-Wigner type interactions, the n-fermion result coincides 
with that of the n-boson result. We may have L m = if we are dealing with fermions with 
some additional degree of freedom such as isospin or colour. Examples of such systems 
include bound states of three quarks or three or four nucleons. 



3 The example of an-almost-harmonic oscillator 

The two step procedure, introduced in the previous section, allows the determination of 
the two-body potential from the knowledge of the spectral information of the n-body 
system. In principle this spectral information for bound states includes both the energy 
levels and the normalisation constants, which represent additional information on the 
wave functions. Both quantities are required to fully determine the potential. Usually the 
energy levels, which are directly related to the mass spectrum of the bound system, can 
be measured with high accuracy. The situation is less satisfactory for the normalisation 
constants, which are not directly accessible to experiment and can only be deduced from 
transition observables. Consequently there is a considerable uncertainty in the normali- 
sation constants. 

It is precisely this situation which we want to study in a first schematic example. 
Given this new n-body inversion procedure, we would like to have an idea of what may be 
learnt about two-body interactions from n-body spectral data. To this end, we perform an 
inversion of n-body spectral information to obtain the two-body interaction in the simplest 
case, where L m = 0. As a basic model then, we consider a bound system of n bosons 
of equal mass m interacting via two-body harmonic oscillator potentials ^. This system 
has only a discrete spectrum and is exactly described by the HCA. The corresponding 
hypercentral potential is also a harmonic oscillator potential, V^ c = (n — 1)~ 1 U, and the 
n-body bound state is obtained by solving a Schrodinger-type equation in the hyperradius 

f d 2 £(£ + 1) p 2 \ . . n 
The energy of this state is and the potential term is simply given by 

^ = S^t£ = £ (26) 



h 2 2 nc hi 



The solution is of the form 




«(A>) = (f) -P(^)^(«.Mf)). (27) 
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where \F\ is the confluent hypergeometric function of the arguments 

a = (2£ + 3 - eb 2 n )/A (21 



and 



/3 = £+~ (29) 



If a = — n, where n is an integer, one obtains a bound state and (27) becomes 



u(p) = N[ £f + ^J^) Ll -i(£)\ ,30) 



where AMs a constant normalising the wave function, and L@ 1 denotes a Laguerre poly- 
nomial. In particular, if a = 0, corresponding to e = e = (2£ + 3)/i£, we obtain the 



ground state w 

w 



It is not our intention to reconstruct the two-body harmonic oscillator potential from 
the corresponding n-body bound-state spectrum because this can be done analytically. 
Rather, we want to simulate the realistic situation taking into account the uncertainties 
in the normalisation constants or energies. Therefore we study the variations of the two- 
body potential due to slight changes in the n-body spectrum. In order to demonstrate the 
essential points in a simple way we restrict ourselves to small modifications of the n-body 
ground states. We alter either the ground state normalisation constant leaving the energy 
unchanged, or the ground state energy, keeping the same ground state normalisation 
constant. All other remaining spectral data of the n-body system are unchanged from 
those of the n-body harmonic oscillator. 

Because of our well defined reference potential it is not necessary to invert this modified 
n-body spectrum via a standard inversion procedure in order to obtain the corresponding 
hypercentral potential. Techniques of supersymmetric quantum mechanics offer a 
more elegant way of achieving the same result in our case. This method is directly related 
to the standard inversion procedure of Gel'fand- Levitan JTTJ and has the additional 



advantage that we can really take into account the whole n-body spectrum. The latter 
point will be very important for the interpretation of our final result because we can be 
sure that all variations found in the two-body potential are generated by the modification 
of the n-body ground-state normalisation constant. There will be no artifacts due to the 
neglect of n-body states, which is unavoidable in a direct application of standard inversion 
methods. 

We alter the n-body spectrum via supersymmetric transformations in two steps, first, 
by removing the harmonic oscillator ground state, and then by inserting a new ground 
state into the n-body spectrum. To remove the ground state from the n-body spectrum, 
we perform the transformation 



d 2 

V=V-2— In M - (32) 
dp 2 



This results in the new potential 



K K p 2 ' ( } 
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The solutions u associated with this new potential may be obtained from equation fl2"T|) 
with the replacements £—>£ + ! and e — > e — thus 



«M=(£f 2 e*p(-!r) I F l(a+ l, /3+ l,(^) 2 ). (34) 



Adding a new state to the spectrum is accomplished by the transformation 

V = V - 2— lnu{l + A / u- 2 {z)dz) (35) 



dp 

where A is a positive constant which controls the bound state normalisation constant of 
the new n-body ground state. 

We need to check that the behaviour of V near the origin does not violate the conditions 
necessary for construction of the corresponding two-body potential. Expanding V — V to 
leading- and next-to-leading order for small p, we find that the change in the potential is 

- 2 < M >" - w +a /; u-w = - * ( |±i - i )( i - 2 ?§±|) +0(A 

(36) 

where ' denotes differentiation with respect to p The p 2 term can immediately be seen 
to cancel the singularity which arises from the removal of the original ground state. 

If we insert a state at the same energy as the ground state of the harmonic oscillator, 
this corresponds to choosing u with a = 0, and the constant in the above expression 
becomes 

1 V- 2^4) = 4^4^^ = ~ (37) 



b 2 n "(3 + l 2 jy 2C + V b 2 n 2C + 5 2C + 1 b\ 

which leaves the potential unchanged at the origin. A change in the energy of the ground 
state will however produce a shift in the potential at the origin. 

If one merely wishes to change the normalisation constant of the ground state one can 
choose u = Uq 1 [I(J and the transformations (p2|) and (p5| ) combine to 

V = V - 2(ln(l + k I" ul(z)dz))\ (38) 
Jo 

where use has been made of the fact that uq is the normalized ground state, and the 
constant K is related to the A of the transformation shown in equation (|35|) by 

-A , , 

(39) 



A+l 

Here the relation between the wave function of the new ground state and that of the 
original ground state near the origin is determined by the choice of k [10| 

limn (p) = y/1 + k limtto(p). (40) 

In our first numerical study we focus on the sensitivity of the hypercentral and cor- 
responding two-body potential to the n-body spectral data at a fixed particle number n 
and orbital angular momentum quantum number / = 0. Using the four-body case as an 
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example, we vary either the size of the normalised four-body ground state wavefunction at 
the origin, or shift the energy by ±5% and ±10%. The changes induced in the potential 
are depicted in Figure 1. While the change in the potential gets larger as the change in 
the energy or normalisation constant increases, it is obvious that even for these relatively 
small changes the potential does not respond linearly. 

We turn our attention to the relative changes in potential that occur if we make the 
same change to the spectra of different number of bodies. In Figures 2(a) and (b) we show 
the change in V£ c induced by a reduction of the bound state normalisation constant by 5 % 
and the lowering of the ground state energy by 5 % without change in the normalisation 
constant respectively for three, four, five and six bodies. The corresponding changes 
induced in the two-body potentials appear in Figures 3(a) and (b). For comparison, we 
include the results in the case of a bound state of two bodies. In Figures 3(c) and (d) 
we show the full new two-body potentials, to illustrate that the change to the oscillator 
potential is not so drastic as to invalidate the use of the HCA. 

Turning first to the changes in V£ c we see that the greater the number of bodies, 
the further from the origin the maximum change in the potential occurs. This is to be 
expected from the form of the HCA: with an increasing number of bodies, the effective 
centrifugal barrier as determined by C increases, and sensitivity to the potential in this 
region is reduced. The change in V£ c decreases as n increases, which provides hope for 
the possibility of obtaining a reliable hypercentral potential from the inversion of n-body 
spectral data. However, this decrease is not transmitted to the underlying two-body 
interaction, as becomes clear on scrutiny of Figure 2. Here the change induced in the two- 
body potential is seen to increase rapidly with increasing particle number, even though 
the change in the corresponding hypercentral potential is decreasing. 

Also worth noting, is that the number of oscillations in the two-body potential in- 
creases with increasing number of bodies. From the form of the transformation fllO|), it 
is clear that p maxima and minima in the hypercentral potential translate into p + k + 1 
extrema in the two-body potential. This has obvious implications for a potential ob- 
tained through inversion, which may have spurious numerically induced oscillations due 
to the neglect of part of the original n-body spectrum. The results clearly indicate the 
limitations on any n-body inversion. 

4 Conclusions 

We have derived a relation for the exact inversion of the underlying two-body potential 
from a given n-body hypercentral potential. In principle this implies that in physical 
systems which are adequately described by the hypercentral approximation, we have ob- 
tained an exact solution of the n-body inverse spectral problem in terms of two-body 
forces. 

We have studied the possibility of extracting information on the two-body potential 
from n-body spectral data in the schematic example of an-almost-harmonic oscillator. It 
turns out that with increasing number of particles the n-body state becomes less sensitive 
to details of the underlying two-body interaction. In other words a small uncertainty in 
the n-body spectral data translates into a dramatic change in the two-body potential. 
This is especially true of any uncertainty in the bound state normalisation constant. 

The harmonic oscillator is a fairly good approximation to some realistic systems, e.g. 
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basic features of light nuclei can be obtained within oscillator models. We are therefore 
justified in transposing our results to realistic systems. In particular they imply that the 
accuracy of the bound state measurements for larger systems must increase considerably 
if we want to get comparable information on the two-body potential. This situation also 
corresponds to an intuitive understanding of the increasing difficulty of extracting the 
contribution to the interaction of a pair of particles, from the total interaction of the 
" ( -" 2 ~ 1 - > pairs, which has been accounted for only in an average way. From this point of 
view it appears to be a further indication of the fundamental conjecture that entropy 
must increase in any many-to-one mapping which is nowhere formally proven. 

Our numerical examples indicate that the inversion of n-body spectra of large systems 
to determine the underlying two-body potential is subject to increasing loss of information. 
However this does not mean that important information cannot be obtained from large 
systems, in addition to that which may be extracted from smaller systems. 

The method remains at present restricted to cases where the hypercentral approxima- 
tion is valid, that is, for weakly correlated systems of particles. An extension of spectral 
inversion for n-body spectra interacting via forces which induce strong correlations be- 
tween particles still remains a formidable problem for the future. 
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Figure Captions 



1. Figure 1: The change produced in the four-body hypercentral and corresponding 
two-body potential if (a) the bound state normalisation constant is altered, and (b) 
the energy of the ground state is altered by ±5% (indicated by dashed and dotted 
lines respectively) and ±10% (the dash-dotted and dash-double dotted lines). 

2. Figure 2: The changes to the hypercentral potential where (a) the bound state 
normalisation constant of the ground state is decreased by 5 % and (b) the ground 
state energy level is lowered by 5 % for three, four, five and six bodies denoted by 
the dashed, dotted, dash-dotted and dash-double dotted lines respectively. 

3. Figure 3: The changes to the two-body potentials (a) and (b), corresponding to 
those in Figure 2 (a) and (b), as well as the total new two-body potential (c) and 
(d). The solid line indicates the (change in) potential for the two-body case. 
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Figure 1(a) 
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Figure 1(b) 
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Figure 2(a) 



Figure 2(b) 
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